Self- Assembly of Patchy Particles into Diamond Structures through Molecular 

Mimicry 



Zhenli Zhang 1 , Aaron S. Keys 1 , Ting Chen 1 , and Sharon C. Glotzer 1 ' 2 ^ 

1 Department of Chemical Engineering and 2 Department of Materials Science Engineering, 
University of Michigan, Ann Arbor MI, 48109-2136 

Fabrication of diamond structures by self-assembly is a fundamental challenge in making three- 
dimensional photonic crystals. We simulate a system of model hard particles with attractive patches 
and show that they can self-assemble into a diamond structure from an initially disordered state. We 
quantify the extent to which the formation of the diamond structure can be facilitated by "seeding" 
the system with small diamond crystallites or by introducing a rotation interaction to mimic a 
carbon-carbon antibonding interaction. Our results suggest patchy particles may serve as colloidal 
"atoms" and "molecules" for the bottom-up self-assembly of three-dimensional crystals. 



I. INTRODUCTION 



The diamond structure is one of the most desirable 
structures from which to make photonic crystals because 
it provides a three-dimensional, complete photonic band 
gap that allows the crystal to diffract light efficiently pQ. 
The current methods of fabrication of these structures in- 
volve direct drilling[2 and layer-by-layer lithography [3], 
which are both top-down approaches and are usually 
expensive and inefficient. Thus the current fundamen- 
tal challenge to the materials community is how to fab- 
ricate the diamond photonic-band gap structure in an 
economically feasible and controllable way that is capa- 
ble of being scaled up to industrial scales. Bottom-up 
self-assembly is a promising strategy, given the right in- 
teractions among the building blocks. There has been 
substantial interest in assembling 3D crystals from col- 
loidal particles [41-HQ]. However, self- assembling a dia- 
mond crystal using a one-component system of colloids 
has not yet been achieved [TT] . Decorating the surface 
of colloids with attractive "patches" suggests a promis- 
ing approach to assemble more complex structures, in a 
potentially controllable and predictable way due to the 
precise interactions between the patches [12]. Many natu- 
ral and synthetic molecules and particles, such as protein 
capsomers in virus shells and nanoparticles with binding 
ligands[T3], can be viewed as patchy particles which may 
serve as programmable building blocks for tomorrow's 
materials [13]. 

Here we propose a model system of patchy particles 
that, soon, may be possible to fabricate based on recent 
experimental findings and theoretical analysis [I5j fl6] . 
By performing Monte Carlo simulations, we show that 
these model patchy particles are capable of assembling 
from a disordered state into a diamond lattice, and we 
quantify the extent to which seeding the system, or 
adding a rotational interaction, can dramatically facil- 
itate the formation of diamond structured assemblies. 
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FIG. 1: Model of patchy particle in our simulations (right). 
Sticky "patches" are shown in pink. Dotted circle indicates 
fourth sticky patch on far side of particle. Related models 
described in text (left). Nelson model in upper left based on 
similar schematic in ref |15j . 



A. Model and Simulation Method 

The model system we consider is illustrated in Fig- 
ure 1 (right), which shows a spherical particle with four 
circular patches arranged in a tetrahedron on the par- 
ticle surface. A pairwise attraction is assigned between 
patches on different particles to allow the particles to as- 
sociate with each other so as to reproduce tetrahedral 
symmetry in three spatial dimensions. The possibility 
of synthesizing this building block experimentally is sup- 
ported by two recent reports. The first is the theoretical 
analysis by Nelson[15 , who proposed a method to make 
tetravalent colloids in which micron-scale spherical col- 
loids are coated by anistropic nanometer-sized objects. 
The ordering of the objects on the particle surface pro- 
duces four disclination defects that may allow the cre- 
ation of tetravalent colloids with chemical or biomolecule 
linkers such as DNA oligonucleotides anchored at the de- 
fect cores, as shown in Figure 1 (upper left). An alter- 
native approach is provided by tetrahedral clusters as in 
Figure 1 (lower left) formed from micron-sized spheres via 
evaporation-driven self-assembly of colloids [16]. By back- 
filling these "secondary building blocks" , it has been pro- 
posed to make spheres in which each of the four colloids 
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within the cluster protrudes slightly beyond the radius of 
the outer sphere, providing circular "patches" of material 
different from that of the outer sphere. To first approx- 
imation, both types of building blocks can be modeled 
by a sphere with four circular patches arranged tetrahe- 
drally. 

The spherical shape of our model patchy particle allows 
the use of a pair potential between patches expressed 
as a hard sphere, square- well potential modulated by an 
angular term [17 
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Here a is the diameter of the particles and set to be 
unity, and e is the depth of the square well potential. 
The function /(Qi, Qj) = 1 only if particles i and j are 
oriented so that the vector joining their mass centers in- 
tersects an attractive patch on both particles; otherwise 
/(Qi,Qj) = 0. The angle J, shown in Figure 1 (right), 
represents the patch size, and A represents the interac- 
tion range. The two values in combination determine the 
directionality and preciseness of the interaction between 
any two patches. In this study we fix the location of the 
patches and choose 5 = ir/6, A = 1.1a which, for a lfim 
diameter particle corresponds to a patch size of 0.26pm 
and an interaction range of 0.1 pm. We leave the effect 
of size, interaction range and disorder in the location of 
the patches for future investigation. 

We perform Monte Carlo simulations using this model 
in a cubic box of fixed size with periodic boundary con- 
ditions. Since we are interested in the diamond struc- 
ture (density p& — 0.6495 for a perfect, box-spanning, 
diamond crystal), we limit our study to p < p&. All 
simulation runs begin from a disordered state at high 
temperature T and are subsequently cooled to the target 
temperature. We investigate three types of systems: sys- 
tems without small diamond crystal seeds, systems with 
seeds, and systems with a modified potential to induce 
the rotation of particles relative to "bonded" neighbors. 
The details of the potential used for the third system will 
be discussed later. 



B. Results and Discussion 

1. Systems Without Seeds 

For the unseeded system, we examine the effect of two 
different cooling rates on assembly. In the "fast" run, T is 
decreased by AT = 0.001 over cycles of 1 million Monte 
Carlo steps (MCS), where T is defined in units of e/k&. A 



MC step is defined here as one attempted MC move per 
particle. In the "slow" run,T is decreased by AT = 0.001 
over cycles of 10 million MCS. A Monte Carlo move for 
a randomly selected particle is randomly chosen to be 
either a translation or a rotation with equal probabil- 
ity. The maximum values of the translation and rotation 
are variable to keep the acceptance ratio at 0.5. De- 
pending on the density, the maximum translation varies 
from 0.02 to 0.04 particle diameters and the maximum 
rotation ranges from 0.05 to 0.09 rad. Although there 
is no simple and direct connection between the number 
of Monte Carlo steps and the physically relevant time 
scales of the system, the purely local, physical moves of 
the simulation provide a reasonable approximation to lo- 
cal diffusive motion. Eleven or twelve independent runs 
are performed to obtain the statistical results shown for 
each density. For the fast cooling rate, only kinetically 
arrested, disordered structures are obtained. However, 
at the slow cooling rate we observe the formation of di- 
amond structures in some of the runs. The fraction of 
successfully assembled diamond structures is 55% (six 
times in 11 runs) at p = 0.40. A typical structure ob- 
tained at p = 0.40 is shown in Figure 2a. To characterize 
the structure we calculate the pair-correlation function, 
g(r), shown in Figure 2b. The correlation function ex- 
hibits peaks at ~ 1.0a, ~ 1.633a, and ~ 1.954a, which 
are the locations of the first three characteristic peaks of 
a perfect diamond crystal. The bond angle distribution 
function, g(6), is shown in Figure 2c. We observe two 
pronounced peaks at ~ 60° and « 180°, which to- 
gether with g(r) confirms the diamond structure for the 
slow cooling rate. 

We calculate the degree of crystallinity cp c and a bond- 
angle correlation function Cb to further quantitatively 
characterize the crystal structures. We define cp c = N s /N 
where A^ s is the number of "solidlike" particles, which we 
define as particles in a locally diamond-like configuration. 
To identify the solid particles, we calculate a local order 
parameter qsm(i) based on spherical harmonics Ysm(^ij) 
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where q3 is a complex vector. If the number of neigh- 
boring particles with which a given particle i is "bonded" 
exceeds three, and if the dot product of c[3(i) and qs(j) is 
less than 0.65, particle i is identified as a solid-like parti- 
cle. The bond-angle correlation function Cb is expressed 

by 



C b = -<q 3 (0-q3(j)*} 



(5) 



We find <p c = 68.9% and C b = 0.779 for the imperfect 
diamond structure in Figure 2a, as compared to corre- 
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FIG. 2: Simulation results of the system without seeds, (a) 
Diamond structure obtained at slow cooling rate for p = 0.40, 
T — 0.11. (b) Pair correlation function g(r) for disordered 
structures obtained with fast cooling rate and diamond struc- 
tures obtained with slow cooling rate. Note absence of char- 
acteristic third peak for disordered structure, (c) Bond-angle 
distribution g(0) for the diamond and disordered structures. 



sponding ideal values of 100% and 1.0, respectively, for a 
perfect diamond crystal. The calculated values are much 
higher than those of the disordered structures, as shown 
in Figure 3, which includes data from several densities. 

Figure 4a shows the time evolution of (p c and Ch for 
the slowly cooled system in Figure 2a. We see that 
both quantities consistently show a sharp transition near 
T = 0.133, which indicates the occurrence of a structural 
transition. In contrast, no obvious change in the corre- 
sponding curves is found in the systems formed via the 
fast cooling rate. A plot of the "phase boundary" be- 
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FIG. 3: Structural information for all systems. SI corre- 
sponds to system without seeds, S2 corresponds to the seeded 
system, and S3 corresponds to the system with "antibond- 
ing" rotational interactions. (Top) Degree of crystallinity (p c . 
(Bottom) Bond-angle correlation function Cb- 



tween diamond and disordered structures determined by 
this method is presented in Figure 4b. 



2. Seeded Systems 

It is well-known that introducing small crystalline 
seeds of the desired structure into a liquid can acceler- 
ate crystallization in molecular and colloidal systems[19j. 
Therefore, it is of interest to ascertain to what extent this 
approach facilitates the formation of the desired diamond 
structure. To introduce the seeds, we freeze eight of the 
512 particles in the simulation box into a single unit cell 
of a perfect diamond structure and let them remain im- 
mobile throughout the simulation. Prior to cooling, all 
other particles except the seed particles are randomly dis- 
tributed in the simulation box. We investigate the same 
two cooling rates as for the unseeded system. 

With the introduction of the seeds, we observe dia- 
mond structures for the fast cooling rate, which did not 
permit crystallization on the time scales of our simula- 
tions in the absence of a seed (Table 1). A typical struc- 
ture at p = 0.40 is shown in Figure 5a. The fraction of 
times the diamond structure is obtained in [TT] indepen- 
dent runs at a fast cooling rate increases from 0% for the 
systems without seeds to 64% for the systems with seeds. 
The same trend is observed for the slow cooling rate; al- 
though the fraction increases less dramatically from 55% 
to 64%, the facilitation is still evident. Table 1 shows 
that seeding facilitates crystallization at nearly all densi- 
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TABLE I: Probability of Observing Diamond Structure for Various Systems. a 





Fast cooling rate 


Slow cooling rate 




SI S2 S3 


SI S2 


,9=0.10 

p =0.20 
p =0.30 
p =0.40 
p =0.50 


0%(0/12) 25%(3/12) 75%(9/12) 
0%(0/12) 42%(5/12) 100%(11/11) 
0%(0/12) 33%(4/12) 100%(12/12) 
0%(0/12) 64%(7/ll) 100%(12/12) 
0%(0/12) 50%(6/12) 100%(12/12) 


25%(3/12) 58%(7/12) 
25%(3/12) 50%(6/12) 
50%(6/12) 75%(9/12) 
55%(6/ll) 64%(7/ll) 
75%(9/12) 67%(8/12) 



a Sl, S2, and S3 are denned in the caption of Figure 3. The percentages in 
the table correspond to N^/N r x 100%, where is the number of runs in which 
diamond structure is observed and N T is the total number of total runs. and 
Nd are indicated in parentheses. 



u 

-a 


















I "*c - 


4 " T = 0.133-^ 


J - Cb 


2 - 


1 1 1 1 1 1 1 



"0 0.5 1 1.5 2 2.5 

MC Cycles (1(f) 



0.16 



0.15 



TO. 14 



0.13 



0.12 



1 1 1 1 


Disordered 


I 1 I 1 


- Q-OS1, Slow 






h-h S2, Fast 






■ S2, Slow 






"S3, Fast 










Diamond 
i.i. 



0.1 0.2 



0.3 
P 



0.4 



0.5 



FIG. 4: (a) Evolution of degree of crystallinity ip c and bond- 
angle correlation function Cb for the system in Figure 2a. (b) 
Plot of maximum T at which diamond structure forms at a 
given density p for systems SI, S2 and S3 as denned in the 
caption of Figure 3. 



ties and cooling rates. Due to the mismatch between dif- 
ferent crystalline domains around the seed, the degree of 
crystallinity and the bond-angle correlation functions for 
the seeded systems are, however, not significantly higher 
than the slowly cooled systems without seeds. 



3. Systems with Particle- Particle "Bond" Rotation 

If we consider the attractive interaction between a pair 
of patches as a type of "bonding" , then all of the bonds 



in our colloidal diamond structure are rotated by 60° or 
180°, as are the carbon-carbon bonds in the staggered 
conformation of ethane. Inspired by this atomistic fea- 
ture of chemical bonding, we introduce an additional po- 
tential energy uij = AU cos(3#) between pairs of bonded 
patchy particles, where is the relative rotation angle 
and AU « 1/cbT. This additional potential energy term 
induces a relative orientation between particles that fa- 
vors the formation of the diamond structure at low T. 

This rotation preference facilitates formation of the di- 
amond structure substantially as compared to the first 
two methods; this is quantified in Table 1. We find dia- 
mond structures for all [12] independent runs at the fast 
cooling rate and density p = 0.40. A typical structure 
is shown in Figure 5b. Further analysis of cp c and Cb 
(Figure 3) indicates the diamond structures have fewer 
defects than those obtained by the first two approaches. 
For example, cp c = 90.9% and C b = 0.922 for p = 0.40, 
indicating the formation of a nearly perfect diamond 
structure. How might such an interaction be induced 
between patchy particles? The rotation of a chemical 
bond in molecules such as ethane is due to antibond- 
ing interactions between the hydrogen atoms on opposite 
ends of the molecule. Here we offer two possibilities for 
mimicking this interaction with particles. The first ap- 
proach involves the use of charges. Assume the patchy 
particle modeled here represents a tetrahedral cluster of 
four spheres fused permanently together. If the colloidal 
spheres carry a uniform surface charge, the repulsive elec- 
trostatic interaction between spheres in different clusters 
will induce a similar rotation, as shown schematically in 
Figure 5c. For 1/im-diameter colloidal spheres, we esti- 
mate that the surface charge needed on each sphere to 
induce an energy difference between two different rota- 
tion states of lk B T at T = 298K is 0.011231 C 2 /m 2 , 
based on a model for a solution of charged colloids or 
proteins in water proposed by Phillies [20]. Considering 
the additional repulsive electrostatic interaction, an at- 
traction energy e of a patch-patch interaction or "bond" 
of roughly 3.0 x 10 _17 J should be sufficient to induce 
formation of the diamond structure. 

The second approach uses a "patterned" patch. In the 
present work, the patches are considered circular and uni- 



5 




FIG. 5: Structure for the systems with seeds and "anti- 
bonding" interactions, (a) Diamond structure for the seeded 
system for p = 0.40, T = 0.11, cooled from T = 0.[14]5. 
(b) Diamond structure for the system with rotational inter- 
actions, for p — 0.40, T — 0.11. (c) Schematic illustration of 
two rotation states of "bonded" tetrahedrally packed cluster 
of four spheres. Top: low energy state. Bottom: high energy 
state. 



form, but they could be made to be nonuniform like the 
linear bi-patch pattern investigated computationally in 
our previous work, which also induces rotation[12]. Thus 
by combining an attractive interaction between patches 
to provide a "bond" between particles with an additional 
interaction that provides bond rotation, colloidal spheres 
can be made to mimic tetravalent carbon atoms. It 
should be noted that the use of seeds and "antibonding" 
interactions apply to thermal systems as well as to sys- 
tems where cooling rate is not a factor, as in suspensions 
of large colloids. 



II. CONCLUSIONS 



In summary, our simulations predict that patchy par- 
ticles with small attractive patches located on the par- 
ticle surface at the corners of a tetrahedron provide a 
promising route to achieve a diamond structure through 
self-assembly. Introducing a small amount of seeds or 
an "antibonding-like" interaction greatly facilitates the 
formation of the diamond structure. Our simulations 
demonstrate that atomic and molecular interactions can 
guide the design of patchy particles as a new generation 
of "colloidal molecules" [HI [HJ [22] to achieve colloidal 
crystal analogues of atomic and molecular crystals. 
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